Entanglement in the Majumdar- Ghosh model 
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We present an analysis of the entanglement characteristics in the Majumdar-Ghosh 
(MG) or J\-Ji Heisenberg model. For a system consisting of up to 28 spins, there is 
a deviation from the scaling behaviour of the entanglement entropy characterizing 
the unfrustrated Heisenberg chain as soon as J2 > 0.25. This feature can be used 
as an indicator of the dimer phase transition that occurs at J2 = J| ~ 0.2411Ji. 
Additionally, we also consider entanglement at the MG point J2 = 0.5Ji. 
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I. INTRODUCTION 

Entanglement has come to be seen as a prime resource for various quantum information 
processing tasks 0]. From another angle, entanglement describes quantum correlations of 
many body systems which on their part are responsible for various interesting physical phe- 
nomena, e.g. quantum phase transitions. In recent years, entanglement has been fruitfully 
used to give an alternate view on quantum phase transitions, especially in low dimensional 
quantum systems (see 0, HI and references therein). Frustrated quantum systems, have 
however been largely left out of the picture in such discussions so far. 

In this Article, we consider the entanglement in a prototypical 1-D frustrated spin system 
- the Majumdar Ghosh or J1-J2 Heisenberg chain (Hf. The spin- 1/2 Hamiltonian is described 
by 

(n.n.) (n.n.n.) 

where the sums run over nearest neighbour (n.n.) and next nearest neighbour (n.n.n.) 
spins, and are spin- 1/2 operators. In what follows, we set Ji — 1 as the energy scale 
and consider antiferromagnetic J 2 > interaction. The ground state properties of this 
system have been studied so far with the use of many methods like exact diagonalization, 
DMRG, field theoretical approach (see and Q for an overview). It is known that the 
model is critical (i.e.jthe spin-spin correlation function decays algebraically) and gapless for 
J2 < J2 ~ 0.2411 0,111. A gap to triplet excitations opens beyond the quantum critical point 
J| accompanied also by the stabilization of a dimerized phase. Interestingly, for J2 = 0.5 
(known as the MG point), the ground state manifold is exactly known. It is spanned by two 
degenerate nearest neighbour dimer states given by 

\R) = (1, 2)(3, 4) . . . (N — 1, N), \L) = (2, 3)(4, 5) . . . (N, 1) (2) 

where denotes a singlet between spins i and j. 




FIG. 1: The model system with nearest J\ and next nearest neighbour interaction J 2 



In the following, we will demonstrate that the entanglement, as quantified by the en- 
tanglement entropy, exhibits characteristic scaling properties in the critical region. We will 
show that the ground state entanglement scales not only with respect to the subsystem size 
n in an infinite system (this possibility was examined carefully in case of models which can 
be solved exactly Q), but in finite systems for fixed n also with respect to the total size 
N of the system. In the gapped phase, a clear deviation from this scaling emerges which 
can be used to identify the critical point of this system. This is to be contrasted with the 
fact that the entanglement measured by standard measures such as entanglement entropy 
or pair-wise concurrence do not detect the critical point of this model, on their own. 

In the second part of this paper, we consider the entanglement at the MG point. This has 
been considered earlier in Ref. y, where it was shown that the nearest neighbour concurrence 
exhibits a jump at the MG point in finite systems. However this jump disappears exponen- 
tially as the total size of the system increases. We derive a simple formula for the change 
of concurrence as the parameter J 2 is made to pass through the MG point. Furthermore, 
we consider a different quantity, the entanglement of a pair of next-nearest neighbour spins, 
as a potential indicator of the MG point. This quantity is maximal at the MG point for 
any size of the system. Furthermore, we analytically show that in the thermodynamic limit 
this entanglement entropy is invariant in the MG manifold, and thus is a robust maximum. 
Interestingly as a byproduct, in the small size limit, the two translationally invariant (or 
'quasi- momentum') states turn out to have entanglement entropies that are local minima in 
this manifold. 



II. ENTANGLEMENT ENTROPY SCALING 



Given a pure state of a quantum many body system, the entanglement of a given sub- 
system with its complement is conveniently measured by the entanglement entropy. The 
entanglement entropy is defined as the von Neumann entropy of the chosen subsystem, i.e. 
if the (reduced) density matrix of the subsystem is p n and Aj are its eigenvalues, then the 
entanglement entropy is defined as 

s (pn) = -^2 Ajlog 2 Xj. (3) 

i 

The entanglement entropy of a block of contiguous spins has been shown to scale differ- 
ently with the block size n, in critical and non-critical ID systems 0. In the thermodynamic 
limit (total system size N — > 00), the entanglement entropy of a non-critical system tends 



2 



to saturate, while it displays universal scaling behaviour for critical systems 

S(n) = c + -log 2 ra, (4) 

where c is a universal scaling constant and c is model dependent. These two constants take 
the values 1 and ir for the isotropic Heisenberg antiferromagnetic chain. The extension of 
Eq.(jlJ) to finite critical systems was given in |10j : 

C i N Tl \ 

S(n,N) = c + -log 2 (-sin(7r-)J. (5) 
For small values of n/N, this equation becomes 

S(n,N) = c 1 + - log2 n -—„*(-) -^^(jf) (6) 

In the rest of this section, we consider the scaling of the block entanglement entropy in 
the ground state of the MG model. The ground state is calculated via exact diagonalization 
(Lanczos method) of systems of up to 28 spins on imposing periodic boundary conditions. 



A. Case 1: Fixed N = 28 

First, let us focus on the scaling of the entanglement entropy w.r.t. the block size n of 
consecutive spins for a system with a fixed number N of spins. This case, for N < 20 and 
J 2 = 0, has been analyzed previously Q. Our results for N = 28 are shown in Fig. [21 
As expected the numerical data for J 2 = (open circles) are well described by Eqs. (jSJ) or 
© describing the saturation of von Neumann entropy for finite N. But there is also good 
agreement of the calculated entanglement for finite J 2 up to J 2 ~ 0.25 with the line given 
by Eq. ©. Fitting the values of S(n,28) to Eq. © yields c = 3.131 and c = 1.017. For 
comparison, the logarithmic divergence of the entanglement in the thermodynamic limit is 
also drawn. 

The correction to this scaling due to the frustrating J 2 in the finite system of 28 spins 
can be most clearly seen for small n. For larger n the finite-size correction (see Eq. and 
the frustration effect show opposite tendencies and cancel each other partially. However, for 
strong frustration near the Majumdar-Gosh point, i.e. for J 2 = 0.49, the frustration effect is 
clearly visible for all n considered. It is thus reasonable to argue that the presence of J 2 7^ 
will produce a saturation of the entanglement vs. n. 

In order to quantitatively characterize the deviation of the entanglement from the scaling 
relation JHJ) at finite frustrating J 2 > 0, the value of \ 2 defined in the standard way, i.e., as 
a square root of the sum of squares of residuals, is calculated and plotted in Fig. El This 
quantity measures the differences between the calculated entropies for J 2 > and the line 
given by Eq. (JSJ). Notably, one can observe a significant deviation from the critical scaling 
for J 2 > J 2 ~ 0.2411. The flat minimum for J 2 w 0.10 is a result of the competition between 
a finite size correction and the correction resulting from the presence of interaction J 2 . 



B. Case 2: Fixed n 

The question arises as to whether the quantum phase transition in the MG model can be 
identified directly from the dependence of "sufficiently local" entanglement measures on the 
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FIG. 2: The scaling of the entanglement entropy of the MG model. The system contains 28 
spins and the length of the subsystem changes from n = 2 to n = 10. The full line represents the 
logarithmic scaling of an infinite system — Eq. Q . The dashed line shows the finite size correction 
given by Eq. ©. Note that the correction resulting due to J2 > is most pronounced for small n. 
For larger subsystem sizes the finite size correction dominates. 



control parameter J 2 , as in certain other models. Indeed, as shown by Osterloh et al. 
the phase transitions in XY spin-1/2 models can be identified explicitly by the dependence 
of the entanglement between two nearest neighbour spins, as measured by the concurrence, 
on the control parameter. On the other hand, single and two-spin entanglement entropy of 
XY models, have been shown to detect quantum phase transitions (see 



The concurrence does not visually detect the dimer phase transition of the MG model, as 
seen m Ref. 0. The single spin entropy is always equal to 1, since the ground state remains 
rotationally invariant (see Section ITTT| . We have checked that entanglement entropies of 
variously chosen subsystems, such as block entanglement of consecutive spins, or blocks of 
next nearest neighbour spins also do not identify the phase transition in this model. The 
reason for this may be that: either local entanglement measures do not detect this phase 
transition or the small system sizes considered do not capture this behaviour [isl |. 

The finite size problem can be circumvented by the usual finite-size scaling techniques. 
Consider e.g. the entanglement entropy of two nearest neighbour spins (n = 2) for different 
total system sizes N. For the isotropic Heisenberg antiferromagnet (J2 = 0) this quantity 
scales as iV -2 as shown in Fig. 0] (the scaling is given by a perfect line with correlation better 
than 0.9999 for N ranging from 16 to 28). This fully agrees with the dominating behaviour 
in Eq.(JHJ). The scaling of blocks of 4, 6 spins are also presented in Fig. 0] 3 • 

Let us focus on the dependence of the scaling of n = 2 entanglement entropy on the control 
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FIG. 3: The least-square deviation ^fx 2 of the calculated values of the entanglement form the 
line described by Eq. [5] (dashed line in Fig. |2J) in dependence on the frustration J 2 . The system 
contains 28 spins and x 2 was calculated for n < 6. Note that \/~x* starts to increase as beyond the 
quantum critical point J| ~ 0.2411. 

parameter J 2 . The correction to the N~ 2 scaling can again be characterized by the value 
of x 2 , describing the difference between the calculated entropies for J 2 > and the straight 
line for J 2 = 0, see Fig. EJ Significantly, there is a clear deviation from critical scaling for 
J 2 > J 2- The phase transition may thus be detected and located by the dependence of ^fx 2 
on the control parameter. 

While the above results show that there is a change in scaling of local entanglement 
entropy around the critical point, the general question posed in this section remains open 
for future discussion. 

III. THE MAJUMDAR-GHOSH POINT 

We now turn to the entanglement at the MG point J 2 = 1/2. In particular, we are going 
to consider two measures of entanglement viz. the concurrence of two spins (which was 
recently analyzed in jij]) and the entanglement entropy of two spins. 

To fix notions, recall that the model Hamiltonian (Eq. ([T|) ) possesses rotational (SU(2)) 
symmetry. It is widely believed that the ground state \g) of 1-dimensional Heisenberg 
antiferromagnets also exhibits this symmetry, i.e. is a total singlet 5 = state. This 
implies that any subset of spins chosen from the whole system is also rotationally invariant. 
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FIG. 4: Finite-size scaling of the von Neumann entropy S(n, N) of the unfrustrated Heisenberg 
chain, i.e. for J2 = 0, for fixed subsystem size n. The system size N changes from 16 to 28 spins. 



This follows from the following simple identity for the reduced state of n arbitrary spins: 
p n = Tr'p g = Tr'U® N p g (Uy N = U® n (Tr'U® N - n Pg (tf)® N - n ) (U^ m = U® n p n (U ] )® n , (7) 

which holds provided p g = \g)(g\ is rotationally invariant, i.e. p g = U® N p g (W)® N for an 
arbitrary single spin unitary operator U (the symbol Tr' denotes the partial trace over the 
unwanted spins). In particular, the reduced state of each individual spin in the ground state 
of the MG model is maximally mixed p\ = 1/2 implying that a single spin is maximally 
entangled with either all or some of the remaining spins in the lattice. Simil arly, any state 
of two spins belongs to the one parameter family of so-called Werner states [15j, that can 
be represented as 

p 2 =p\^)(y-\ + L-lt, (8) 

where denotes a singlet and p = — (4/3) (S^ • Sj) is just the (rescaled) isotropic corre- 
lation function of the involved spins. 

The entanglement between two spins, characterized by the concurrence C(p2), can 
be checked to be given by a simple formula for Werner states 

C(p 2 ) = max (0, |p - i) . (9) 

Thus, two spins are entangled with each other as long as the correlations between them are 
sufficiently antiferromagnetic, (Si ■ Sj) < —1/4 [20]. The entanglement entropy (see Eq. 
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FIG. 5: The deviation (measured by y x 2 ) from the straight line scaling for J2 
upper line in Fig. |3J for bigger values J2. 



and n = 2 (the 



©) of two spins in a Werner state, on the other hand, is given by the relation 
S(p 2 ) = 2 - log(l + 3p) - 3^ log(l - p). 



(10) 



Returning now to the MG point, one of the states \R) or \L) is realized as the ground 
state, with broken translation symmetry, in the thermodynamic limit. In these states, the 
nearest neighbour concurrence is either or 1 depending on whether the considered pair 
resides on the same or different singlets. The average nearest neighbour concurrence in both 
these states is C(\R(L))) av = 1/2. For finite chains however, in the absence of an additional 
symmetry breaking field, it is more natural to revert to an orthogonal "qubit" basis of the 
ground state manifold, which can be chosen to be the eigenstates of the momentum operator 



l±> 



/a 



--(\R)±\l)), 



where the normalizing factors are 



n ± = 2(l±x), x = (R\L) = (-l) N ' 2 2 1 - N ' 2 



(11) 



(12) 



and x is the overlap of the two dimer states. On traversing the MG point from left to right, 
the ground state changes from |+)(|— )) to |— )(|+)) for a translationally invariant system 
with even (odd) N/2 (this is related to Marshall's sign law 0]). It is thus interesting to 
characterize the entanglement in the momentum basis |±). 
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FIG. 6: Entropy of a pair of n.n.n. spins in the ground state manifold at the MG point 



First, consider the entanglement of two nearest neighbour spins, say 1 and 2 (the same end 
result holds for all nearest neighbour spins in these states, due to translational invariance). 
The parameter p could be determined by calculating the corresponding correlation functions. 
Equivalently, we calculate the form of the state pnp) of these two spins explicitly. Notice 
that the two spins are bound into a singlet in the state \R), while they belong to different 
singlets in the state \L) thus yielding a maximally mixed reduced state. Hence 



= Tr'|±><±| = —Tr'(\R)(R\ ± \R)(L\ ± \L)(R\ + \L)(L\) 



1 



n 



— ((i±2x) |*-)(*- fl/4). 



(13) 

The entanglement between two neighbouring spins in the symmetric and antisymmetric 
states is determined by the parameters 

1 + 2r 

(14) 



P± = 



2(1 ±x) 



For N > 6, both states are entangled (p± > 1/3), and so the difference in the concurrence 
between them can be calculated from Eq.Q 



AC 



cipiV: 



C(p 



12 , 



:(P+~P- 



3x 



2(1 



x z 



(15) 



The absolute value of this expression gives the "jump" in the nearest neighbour concurrence 
on traversing the MG point from left to right (the sign of this difference depends on N, which 



S 



is consistent with the ground states in the vicinity of the MG point.). This quantity has been 
proposed to be an indicator of the MG point in the concurrence diagram in Ref. 0- However, 
for large N, AC approaches zero exponentially. Additionally, two spins that are not nearest 
neighbours are not entangled with each other, since the correlation function drops rapidly 
with distance. Thus a more general measure considered in jij], viz. the total concurrence 
being the sum of concurrences of all pairs contains only one non-zero contribution coming 
from nearest neighbours and hence cannot detect the MG point. 

The question thus arises concerning other indicators of this special point. For the infinite 
system, a simple candidate could be the dimer order parameter given by the difference 



1 

N 



(16) 



At the MG point, d takes the value 3/8. Considering the notion of dimerization, one 
could naively assume that this is the largest possible value as it corresponds to exact dimers. 
However, the states \R), \L) are not eigenstates of the dimer operator and as such d cannot 
take on extremal values for these states. Physically, one could expect that high but not 
perfect antiferromagnetic correlations on one bond supplemented by slightly ferromagnetic 
correlations on the other could create a larger value of the dimerization. Using DMRG 
techniques, it has been shown that indeed the maximum dimerization does not occur at the 
MG point but at J 2 « 0.5781 ja). Additionally, this quantity depends largely on the working 
basis: for translationally invariant states like |±), the value of this parameter is always zero. 
Thus the dimer order parameter also does not distinguish the MG point satisfactorily. 

The entanglement entropy of a pair of next nearest neighbour spins is a much more 
appropriate quantity that distinguishes the MG point. Due to degeneracy at the MG point, 
the general form of the ground state is 



I*,) = ■ = (cos e -\R) + sm 9 -\L)). (17) 

V 1 + xsm v \ £ £ ' 



The parameter for a pair of n.n.n spins in the above state is given by 

x sin 6 

Pi,i+2 = -7— r-j. 18 

1 + x sm d 

For finite systems, the dimer states \R), \L) obviously maximize the entropy of entangle- 
ment which is equal to 2. Moreover, as the size of the system increases, the dependence of 
the entropy on 9 flattens out to the maximal possible value exponentially fast, since x — > 
(see Fig. EJ). This further justifies the choice of this quantity as a universal indicator of the 
MG point. Interestingly, the momentum states |±) are distinguished as local minima in the 
ground state entanglement diagram. In the wider range of values of J2, the pair entropy 
of the MG point is indeed uniquely distinguished (see Fig. for the n.n.n. entanglement 
entropy for 16 spins) as the sole maximum in the diagram. 



IV. CONCLUDING REMARKS 

We have focussed on the entanglement properties in the Majumdar-Ghosh model. Based 
on data from numerical calculations of finite chains (up to 28 spins), we have discussed the 
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FIG. 7: The n.n.n. pair entanglement entropy for N = 16 spins. The two jumps which are finite 
size effects, connected to energy level crossings for this system. These discontinuities vanish as the 
system size grows. 

scaling of the entanglement entropy of blocks of consecutive spins and shown that it can 
be used as a tool to identify the quantum critical point of this model. In contrast with 
other numerically studied models, the critical behaviour of the system does not manifest 
itself directly in the dependence of "local" entanglement on the control parameter, for the 
considered system sizes. However, the transition from the critical gapless phase to the non- 
critical gapped phase appears in the characteristic change in scaling of the local entanglement 
measures w.r.t. total system sizes. Furthermore, we have shown that the Majumdar-Ghosh 
point of this model can be identified as a maximum in the dependence of next nearest 
neighbour pair entanglement on the control parameter. 

In the end, we would like to add that one can heuristically consider the entanglement 
entropy of the lower rail of spins with the upper rail of the considered system (Fig. as 
a natural candidate for distinguishing the phases of this model. Indeed, the dimer phase 
is expected to be characterized by enhanced correlations between the lower and upper rails 
of spins. Once again however, for the system sizes considered there is no "characteristic 
change" in the dependence of this quantity on the control parameter. The results have not 
been provided in this Article, since they resemble the results presented in Fig[7| At the MG 
point, again the "momentum" eigenstates reside in local minima, as in Fig. |U1 
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The systems sizes studied in Ref. Q range from A = 6 to 12, while we have densely checked 
for up to N = 18. 

For these and larger blocks, effects of terms of the order O (n/N) become important. 
Interestingly, {Si ■ Sj) = — 1/4 is the classically optimal antiferromagnetic correlation corre- 
sponding to the Neel state. So, enhanced antiferromagnetic correlations are a manifestation 
of quantum entanglement. 
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